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We calculate next-to-leading-order (NLO) correction to the pion electromagnetic form factor 
at leading twist in the fcr factorization theorem. Partons off-shell by are considered in both 
quark diagrams and effective diagrams for the transverse-momentum-dependent (TMD) pion wave 
function. The light-cone singularities in the TMD pion wave function are regularized by rotating 
the Wilson lines away from the light cone. The soft divergences from gluon exchanges among initial- 
and final-state partons cancel exactly. We derive the infrared-finite fcT-dependent NLO hard kernel 
for the pion electromagnetic form factor by taking the difference of the above two sets of diagrams. 
Varying the renormalization and factorization scales, we find that the NLO correction is smaller, 
when both the scales are set to the invariant masses of internal particles: it becomes lower than 40% 
of the leading-order (LO) contribution for momentum transfer squared > 7 GeV^. It is observed 
that the NLO leading-twist correction does not play an essential role in explaining the experimental 
data, but the LO higher-twist contribution does. 

PACS numbers: 12.38.Bx, 12.38.Cy, 12.39.St 



I. INTRODUCTION 

The kx factorization theorem [ll-Q in perturbative QCD (PQCD) has caught a lot of attention recently (see 0, ll] 
for its recent applications). How to derive a hard kernel beyond leading order (LO) in the kx factorization theorem 
is one of the important subjects. We have proposed the prescription that partons in both QCD quark diagrams and 
effective diagrams for transverse-momentum-dependent (TMD) hadron wave functions are off mass shcU by fc| d. 
Since the gauge dependence cancels between the two sets of diagrams, their difference defines a gauge-invariant hard 
kernel. The gauge invariance of the kx factorization theorem has been investigated in [lol and proved in p^ . 
The light-cone singularities, appearing in the naive definition of TMD wave functions, are regularized by rotating the 
Wilson lines away from the light cone [l3l - [T5| . Following the above prescription, the next-to-leading-order (NLO) 
correction to the pion transition form factor associated with the process nj* — >■ 7 has been calculated at leading twist 
prj . Both the large double logarithms as hi^ kx and as In^ x, x being a parton momentum fraction, were identified. 
The former is absorbed into the pion wave function and summed to all orders in as by the kx resummation, and the 
latter is absorbed into a jet function and summed to all orders in as by the threshold resummation. 

In this paper we shall extend the framework to the more complicated pion electromagnetic form factor associated 
with the process ttj* — )- tt. Our work represents the first NLO study of this quantity at leading twist in the fcy 
factorization theorem. The NLO correction to the pion form factor has been evaluated in the coUinear factorization 
theorem by several groups [l6l - [23 |. For 777* — >■ tt, not only the coUinear divergences from gluon emissions collimated 
to the initial- and final-state pions, but also the soft divergences from gluon exchanges between the two pions exist. 
Therefore, it is nontrivial to demonstrate that the coUinear divergences in the quark diagrams are cancelled by those in 
the pion wave functions, and the soft divergences cancel among themselves. The computation of the box and pentagon 
diagrams for the pion form factor is also technically involved. The renormalization scale /i and the factorization scale 
/if are introduced by higher-order corrections to the quark diagrams and to the effective diagrams, respectively. We 
shall show that the NLO correction to the pion form factor is smaller, when both the scales are set to the invariant 
masses of internal particles as postulated in 0, [2^ : it becomes lower than 40% of the LO contribution for momentum 
transfer squared > 7 GeV^ . The variation of the pion form factor with /.j and {if gives an estimate of the theoretical 
uncertainty in our analysis. 

Another related subject concerns the shape of the pion distribution amplitudes. It has been pointed out [2^ that 
the recent BaBar data of the pion transition form factor ^^-^(Q^) [13] at low (high) indicate an asymptotic (flat) 
leading- twist, i.e., twist-2 pion distribution amplitude. However, a pion distribution amplitude ought to be universal. 
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These seemingly contradictory observations have been reconciled in the fcy factorization theorem [26j : the increase of 
the measured Q'^Ft^^[Q'^) for Q"^ > 10 GeV^ was explained by convoluting a fcy-dependent hard kernel with a flat pion 
distribution amplitude. The low data were accommodated by including the threshold resummation of as In x, 
which provides a strong suppression at the endpoints of x. The suppression becomes stronger as decreases. That 
is, the twist-2 pion distribution amplitude multiplied by the threshold resummation factor behaves like the asymptotic 
model at low effectively. In this work we shall show that the asymptotic models for both the twist-2 and two-parton 
twist-3 pion distribution amplitudes also lead to results in better agreement with the data of the pion form factor at 
low Q^, compared to those from the non-asymptotic models. It is found that the NLO twist-2 correction does not play 
an essential role in accounting for the data, but the contribution from the LO two-parton twist-3 pion distribution 
amplitudes does, a conclusion the same as drawn in [28|. 

In Sec. II we calculate the 0{a^) QCD quark diagrams, the 0{as) effective diagrams for the TMD pion wave func- 
tion, and their convolution with the 0{as) hard kernel, considering partons off-shell by kjn. Since the kr factorization 
theorem is appropriate for QCD processes dominated by contributions from small x [9|, we shall keep only terms 
in leading power of x. Taking the difference between the two sets of diagrams, we derive the NLO fcy-dependent 
hard kernel for the pion form factor. Section III contains the numerical investigation, in which we obtain the partial 
contributions to the pion form factor from different twists and orders, and the relative importance of the LO and 
NLO twist-2 contributions. We also examine the dependence of our results on the renormalization and factorization 
scales, and on the shape of the pion distribution amplitudes. Section IV is the conclusion. 



II. NLO CORRECTIONS 



In this section we compute the 0{a'^) quark diagrams in QCD and the 0{as) effective diagrams for the pion 
wave function in the Feynman gauge, and then derive the NLO hard kernel for the pion electromagnetic form factor 
in the fc^ factorization theorem. The momentum Pi {P2) of the initial-state (final-state) pion is chosen as Pi = 
{P^,0,Ot) {P2 ~ (0,^2"' Ot))- The anti-quark q carries the momentum ki = (xiPj^, 0, ki-r) in the initial state and 
^2 — {0,X2P2^ ,^2t) in the final state, as indicated by the LO quark diagrams in Fig. [1] The parton virtuality k'j^ 
regularizes the collinear and soft divergences into the logarithm lnfc|,. We shall focus only on Fig. [IJa), since the 
NLO correction to Fig. [IJb) can be obtained from that to Fig. [IJa) by exchanging the kinematic variables of the two 
pions. The NLO corrections to the other two LO quark diagrams with the virtual photon attaching to the anti-quark 
line are then obtained via the variable exchanges between the quark and the anti-quark. According to the hierarchy 
^ xQ^ ^ fc|, in the small-x region, we keep only those terms that do not vanish in the x — ?> and fc^ — > limit. 

Figure [TJa) leads to the amplitude 

H^''>{xi,kiT,x2,k2T,Q ) = ~g Cf—^=-:t 7^ TTWTT r-k ' 

= n^r A7m A) 

^ ^Q^XiX2Q^ + \klT-\^2Tn' ^' 

where Nc is the number of colors, Cp is a color factor, A75/\/2iVc and 75 A/v^^lVc are the twist-2 spin structures of 
the initial- and final-state pions, respectively, and = 2Pi ■ P2 is the momentum transfer squared from the virtual 
photon. To reach the last line in Eq. ([ij, we have applied the hierarchy xiQ^ ^ k^rp to the internal quark propagator. 
The denominator xiX2Q^ + |kiT — k2Tp comes from the virtuality of the LO hard gluon. The |kiT — k2Tp term may 
not be negligible compared to XiX2Q^ in the small- a; region, so it is retained. 
The NLO hard kernel _ff is defined, in the kx factorization theorem, by @ 

H'^'\xukiT,X2,k2T,Q^) = G(l)(.Ti,fciT,a:2,fc2T,Q') 

- J dx[(f k[rj.^^^\xi , kiT] x[, k[rp)H^°^ {x[,k[rp, X2, fc2T, Q^) 

-J d4d2^^^^(°)(a;i,fciT,x^,fc^r,g2)$(i)(a;2,fc2T;4,fc2T), (2) 

where G^^^ denotes the NLO quark diagrams associated with Fig. [TJa), and collects the 0{as) effective diagrams 
for the twist-2 quark-level wave function [2^ 



<!>{xi,k,T;x[,k[r) = J ^^e-"'^^i+'^"+^''--y-(0|g-(2/)VF,(n)t/„,,,oW^oW A-75g(0)k(Pi - ^1)9(^1)). (3) 
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(c) [d) 
FIG. 1: Leading-order quark diagrams for 7r7* — s- tt with ® representing the virtual photon vertex. 



In the above expression y = (0,y tYt) is the coordinate of the anti-quark field g, n_ = (0, 1,0t) is a nuh vector 
along P2, k(-Pi — ki)q{ki)) is the leading Fock state of the pion, and the Wilson line Wy(n) with 7^ is written as 



Wy (n) = P exp 



-ig I dXn ■ A{y + \n) 
'0 



(4) 



The two Wilson lines Wy{n) and WQ{n) are connected by a vertical link /,i;j,,o at infinity [30, HH. Equation ([3]) will 
generate additional light-cone singularities from the region with a loop momentum collinear to n_, as the Wilson line 
direction approaches the light cone, i.e., as n — )■ n_ [l3|. That is, serves as an infrared regulator for the light-cone 
singularities. 



A. NLO Quark Diagrams 

We first calculate the NLO corrections to Fig. [lla) in the fc-r factorization theorem, which come from Figs. [21 [3l 
and [H for the self-energy corrections, the vertex corrections, and the box and pentagon diagrams, respectively. The 
ultraviolet poles are identified in the dimensional reduction f32| in order to avoid the ambiguity from handling the 
matrix 75. To simplify the expressions, we define the following dimensionless ratios 

- Q2 ■ (5) 

Terms suppressed by powers of cc or 5 will be dropped in the NLO corrections. 
The self-energy corrections to the four external quarks in Figs. [Ha)-[21Jd) give 

where 1/e represents the ultraviolet pole, fi is the renormalization scale, and is the Euler constant. The collinear 
divergences in Figs. [llja)-[21[d) are regularized into the infrared logarithms In S. The self-energy correction to the 
internal quark in Fig. ^e) leads to 

^(r) ^ (^i + 1, + 2) ij(0). (8) 
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FIG. 2: Self-energy diagrams. 



FIG. 3: Vertex-correction diagrams. 



The ultraviolet poles in Eqs. ([5]) and ([7]) arc half of that in Eq. ([S]), because an additional factor 1/2 is associated with 
an external particle. The internal quark is off-shell by the invariant mass squared xiQ^, which replaces the arguments 
of the infrared logarithms in Eqs. 1^ and ([7]). The self-energy correction to the hard gluon, summing the quark-, 
gluon- and ghost-loop contributions in Figs, [^rf)-^!), is written as 



'^2f+2g+2h+2i - I 3 ^ 



(9) 



with N f being the number of quark flavors. In this case the logarithm depends on the invariant mass squared Si2Q^ 
of the hard gluon. 

The results from the five diagrams Figs. [3][a)-[3l^e) for the vertex corrections are summarized as 
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The ultraviolet poles in the sell-energy and vertex corrections are summed into 
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for Nc = 3, which is the same as obtained in |22| . This pole term determines the renormalization-group (RG) evolution 
of the coupling constant a^. 

The loop correction Gg^j'' to the photon vertex does not contain an infrared logarithm, because of the sequence of 
the 7-matrices in the fermion flow. The correction Ggj^' to the upper gluon vertex contains only the infrared logarithm 

associated with the final-state pion. The other end of the radiative gluon attaches to the internal quark line, so G^i^^l 
depends on xiQ^, and does not generate a double logarithm. The correction to the lower gluon vertex in Fig. [2Ic) 
depends on the momentum transfer squared 5i2Q^ from the LO hard gluon, and on the two infrared regulators 5i and 
82 from the incoming and outgoing partons. The expression of Gg^^ is symmetric under the exchange of 61 and 82 as 
expected. The ln((5i2/(5i) ln(5i2/^2) term is not a large double logarithm in the small-x region, where h\6i In (52 arises 
from the soft region of the radiative gluon. It is obvious from the convolutions of ^'^^^ and H'^'^^ in Eq. ^ that this 
soft logarithm can not be cancelled by the infrared logarithms in the effective diagrams: there is no chance for In i5i 
associated with the initial-state pion and ln<52 associated with the final-state pion to appear in a product. Therefore, 
cancellation must occur within this type of soft logarithms as demonstrated below. 

The triple-gluon vertex correction Gg^ also depends on the scale 5i2Q^ , and contains the same infrared logarithm 



hi{5'l2/ 8182) as in Gg^'' but with different color factors. Their sum is proportional to Cp'- —^/Nc + Nc = 2Gf, 
implying the factorization of this infrared logarithm in color fiow. It is then possible that it will be cancelled by the 
corresponding logarithm in the effective diagrams, whose color factor is also Cp- This combination of quark diagrams 
for achieving factorization also takes place elsewhere. Another triple-gluon vertex correction in Fig. [3{e) involves the 
invariant masses squared of the virtual quark and of the hard gluon as shown in Eq. (|14p . At small x, the quark q 
in Fig. [3](e) is energetic, leading to the collinear logarithmic enhancement ln(a;i/i52), and the hard gluon has a small 
invariant mass, leading to the soft enhancement \n(xi/8i2). Their overlap then gives rise to the important double 
logarithm ln(a;i/(52) ln(a;i/(5i2) in Ggg''. This double logarithm can be reexpressed as 



21n^ln^ 
82 812 



2 Xl 



8u 



In-Zl+ln-— _ln^^, (16) 

O2 Xl 82 

where the first term, known as the Sudakov logarithm 0, [s^ , is absorbed into the final-state pion wave function. 
The Sudakov effect from resumming this double logarithm suppresses the contribution from the small k2T region, 
i.e., the region with a large impact parameter The second term \n^{8i2/ xi) « In^ X2 also exists in the collinear 
factorization theorem with kx being integrated out f3^-'36'| . This threshold logarithm is significant at small X2 , where 
the hard gluon approaches mass shell [37.] . The third term in Eq. ()16|) is less important in the small-x region. 

There are four box diagrams and two pentagon diagrams as displayed in Fig. SI Note that Figs. ID^c), Hfe) and 
mjf) also contain the corrections to the LO quark diagram Fig. [TJb). As explained in [ll|, the two-particle reducible 
diagrams, Figs.^a) andlUe), give power-suppressed contributions at small x. The reason is simply that a collinear or 
ultraviolet loop momentum increases the virtuality of the LO hard gluons. A soft loop momentum does not change the 
power-law behavior of the hard kernel, which is, however, not a leading region due to the soft cancellation mentioned 
before. The expressions for the diagrams other than Figs.SlJa) andSl^e) are collected below: 
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with being the LO hard kernel from Fig. Wlh) 



fV(0)/ , r)2\ 2^ /Pi /P27^ /P2 

Q'^{xiX2Q^ + kiT - k2Tr) 



(21) 



The above expressions are free of ultraviolet divergences. The correction G^^^ involves only ln(5i, if focusing on the 
correction to Fig. [Ija), and the important double logarithm ln(5ilnxi: Ini^i denotes the coUinear enhancement, and 
Inxi denotes the soft enhancement in the small xi region. Similarly, the above double logarithm can be reexpressed 



2 In ^1 In = In^ 5i + In^ Xi — It? , 



(22) 



for the separation the Sudakov logarithm and the threshold logarithm. Figures. |4l[b), lljd), andlH^f), with the soft 
radiative gluons attaching to the incoming and outgoing partons, generate the soft logarithm In (5i In (52. These soft 
logarithms cancel among Figs. [Sjc), IH^b), |4l[d), andSJ^f), for the reason that soft gluons do not interact with color- 
singlet objects like pions. 

The sum over all the NLO quark diagrams associated with Fig. [IJa) gives 
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for Nf — 6. The NLO correction to Fig. [TJb) can be obtained from the above expression by performing the exchanges 
Xi O X2 and Si ■<r^ 62- 



B. NLO Effective Diagrams 

We compute the convolutions of the NLO wave functions with the LO hard kernel over the integration 
variables x' and kip, 

= J dx[d^k[T<i>['HxukiT;x[,k[p)H'^''Hx[,k[T,X2,k2T,Q^), 
i?^"^^*!" = J dx'2d^k'2TH^''\xi,kiT,x'2,k'2T,Q^)^i'\x2,k2T;x'2,k'2T), (24) 




with a choice of 7^ for the Wilson hne direction to regularize the light-cone singularities. The NLO wave 
functions then depend on through the scale = 4(n • Pi)^/|n^| or Cl = 4(n • P2)^/|?t-^|, which is regarded as a 
factorization-scheme dependence. This dependence, entering a hard kernel when taking the difference between the 
quark diagrams and the effective diagrams, can be minimized by adhering to a fixed n'^. Note that the soft subtraction 
factor introduced in [l^ is not necessary here, which is needed for a choice of = 0. 
The self-energy corrections from Figs.lSfa) and[SlJb) are written as 
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whose expressions are similar to those from the quark diagrams but with the factorization scale /if. The contribution 
from the box diagram in Fig. [5jc) is power-suppressed in the small x region: 



$(1) ® = ^ $(1) = 0. 



(27) 



The sign of the plus component n+ of the vector n is arbitrary, which could be positive or negative (n has a 
positive sign, the same as of P2 )• Choosing n+ < 0, i.e., < as in SHUi, Fig.IHd) leads, in the small x region, 
to 
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which reproduces the Sudakov logarithm in the form of In {(^ /k'^). As computing the convolution of "I>i ' with H'^^\ 
the momentum fraction appearing in the hard kernel should be restricted between and 1. The expression for Fig.[5][e) 
is given, in the small x region, by 
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where terms vanishing with /c^ have been dropped. It is observed that Fig.[5je) also generates a double logarithm, 
whose importance is attenuated by the small x. 
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The result from Fig. ^{) is similar to that from Fig. ^d), but with the replacement of P — fc by fc, i.e., ( by x(. 
Keeping terms which do not vanish with /c|, — >■ 0, we have 
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where the double logarithm is further attenuated by x^. It should disappear, after combined with the contribution 
from Fig.lSJg), since such a double logarithm is absent in the NLO quark diagrams. The same variable transformation 
relating ^y-^"* to ^^^^ is not applicable to ^g^"*, for the latter involves the nontrivial convolution with H^^\ Retaining 
terms which are finite as kx — >■ 0, Fig. [Sfg) leads to 
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It is easy to see the cancellation of the double logarithms between Eqs. (150]) and (PT|) . 
At last, we include the self-energy corrections to the Wilson lines, namely. Fig. [5jh-j) 
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Summing all the above 0{as) contributions, we obtain 
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The resultant anomalous dimension of the pion wave function is the same as derived in the axial gauge (for a recent 
reference, see [Slj). 



C. NLO Hard Kernel 



We derive the NLO hard kernel for the pion electromagnetic form factor in the fcr factorization theorem by taking 
the difference between the quark and effective diagrams. Note that ag appearing in Eqs. (|23|) and ([33]) denotes the 
bare coupling constant, which can be rewritten as 

as = asim) + SZ(pi)as{pi), (34) 

with the counterterm SZ being defined in the modified minimal subtraction scheme. We insert Eq. ([M)) into the 
expressions of the LO and NLO quark diagrams, and of the NLO effective diagrams. The LO hard kernel H'-°^ 
multiplied by SZ then regularizes the ultraviolet pole in Eq. (j23p . The ultraviolet pole in Eq. (1331) is regularized by 
the counterterm of the quark field and by an additive counterterm in the modified minimal subtraction scheme. For the 
infrared logarithms, there exists a cancellation between the two-particle reducible quark diagrams and the two-particle 
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reducible effective diagrams. That is, 0^2! 2b ~ ^^'^^ ® *^a^6 ^2^c2d ~ *^a^fc ® 7?^°-* are infrared finite. The self-energy 
diagrams for the internal quarks and gluons are free of infrared divergences. The collinear logarithm represented by 
ln(5i (ln(52) cancels between the two-particle irreducible diagrams and the convolution ^^^g\®U^'^^ 

The infrared-finite /c^-dependent NLO hard kernel for Fig. [IJa) is given by 

ti^ ' = -T 7=rT ~ — -In xi -f — Inxi lnx2 - in 0\2 + —r mxi lndi2 

47r \ 4 y-^ 4 8 4 

13, 31 23, ^ 1, ^ 57r2 53\ .q) 

- y In a:i + — In X2 + y In (5i2 + - In 2 + — + — j . (35) 

We have chosen |n+| — , which renders ("^ = Q^, to avoid creating an additional large logarithm like \n{(^^/Q^). 
The Sudakov logarithm \n^{Q^/k^) in Eq. ([^^ has been removed by that in the effective diagrams, but the threshold 
logarithms In remains in which can be absorbed into a jet function [s^]. Because there is no end-point 

singularity in the pion form factor at twist-2, we shall not perform the factorization of this threshold logarithm 
here. However, an end-point singularity does exist in the contribution from the two-parton twist-3 pion distribution 
amplitudes, so the factorization and the threshold resummation for the jet fimction needs to be performed. 

We have checked that Eq. ([55)1 is gauge-invariant by repeating the NLO calculations in an arbitrary covariant gauge. 
Applying the Ward identity to the gauge-dependent terms [lll[l2|> it is easy to show that the gauge dependence cancels 
between the quark diagrams and the effective diagrams. The NLO hard kernel H'-^^ for Fig. Hfb) can be obtained 
from Eq. ([35]) by exchanging xi and X2, and by substituting i/'"-* in Eq. ((2T|) for H'^^\ The NLO hard kernels for 
the other two LO quark diagrams with the virtual photon attaching to the anti-quark line can be obtained from H^^^ 
and H^^^ by exchanging x and 1 — x. 



III. NUMERICAL ANALYSIS 



In this section we evaluate the NLO correction to the pion electromagnetic form factor in the kr factorization theo- 
rem numerically. Since the Sudakov suppression was derived in the impact parameter space, we write the factorization 
formula as a convolution in the impact parameters 61 and 62 conjugate to kir and k2T, respectively. The Sudakov 
factor, resulting from the summation of the double logarithm as In^ kx to all orders, describes the perturbative kx 
dependence of a pion wave function. We shall not consider the intrinsic kx dependence of a pion wave function pro- 
posed in [13] here. To simplify the expression of the NLO hard kernel in the &1-62 space, we adopt the approximation 
ln(5i2 ~ ln(xia;2)- This approximation does not modify the power-law behavior of end-point contributions from small 
X, and affects numerical outcomes only slightly. We first test the following asymptotic models for the twist-2 pion 
distribution amplitude (j)^ and the two-parton twist-3 pion distribution amplitudes 0^'"^, 

<j>^ix) = -%=x(l-x). 



(l>i{x) = -^(l-2x), (36) 



where the pion decay constant is taken as /^r = 130 MeV. 

The first issue concerns the choice of the renormalization scale fj, and the factorization scale /Xf in order to minimize 
the NLO correction to the pion form factor. The results are summarized in Fig. [6l which displays the ratio of the 
NLO contribution over the LO one as a function of Q^. For the first choice, /j,f is set to the hard scale t defined in 
the PQCD approach to exclusive processes based on the kx factorization theorem 0, [25j 

t = max(V^g, 1/61, 1/62)- (37) 

Then we utilize the freedom of choosing fi to diminish all the single-logarithmic and constant terms in the NLO hard 
kernel, which is found to be 

, / N /53 57r2\ 5/42 11/24 /wV^^ /qQ^ 

t,{fif)^exp[— + —]x^' X2' — Mf. (38) 



42 504 y ' ^ \Q 



10 



1.0 




0.3- 



0.2 J—, ■ 1 ■ 1 : 1 : 1 : T" 

20 40 60 80 100 

Q' (GeV") 



FIG. 6: Ratio of the NLO correction over the LO contribution to the pion form factor with the scale ts(^f) being defined in 
Eq. (|38j- 
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FIG. 7: Pion form factor from different orders and twists for the asymptotic pion distribution amphtudes in Eq. (|36|) . The 
experimental data are taken from (4L] and (42]. 



This scale choice gives a larger NLO correction: the NLO correction becomes less than 40% of the LO contribution 
only when is higher than 55 GeV^. The second choice corresponds to /if = Q, which is greater than t, and 
H = tg, which removes all the single-logarithmic and constant terms in the NLO hard kernel. The convergence of the 
NLO correction is improved: the NLO correction becomes less than 40% as > 17 GeV^. The third choice is the 
conventional one adopted in the PQCD approach to exclusive processes, with both fi and /if being set to t. It turns 
out that this simple choice works better: the NLO correction becomes less than 40% as > 7 GeV^. The above 
analysis confirms the prescription postulated in [HI that both the renormalization and factorization scales are set to 
the invariant masses of internal particles. 

After fixing the scales, we calculate the LO twist-2, NLO twist-2, and LO two-parton twist-3 contributions to the 
pion form factor, and present the results in Fig. [T] Note that the last piece of contribution involves two two-parton 
twist-3 pion distribution amplitudes, so it is suppressed by a power of l/Q^- Figure [7] indicates that the LO twist-2 
contribution is quite small compared to the experimental data in the whole considered range of . The NLO twist-2 
contribution is also small, and helps to explain the data only at ~ 1 GeV^, where the perturbation theory may 
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FIG. 8: Contributions to the pion form factor from different orders and twists for the non- asymptotic models in Eq. (|39p . 



not be reliable: the quick growth of our results at very low is attributed to the presence of the infrared Landau 
pole in the strong coupling constant. To study the pion form factor in this region, a soft contribution is needed, and 
an analytic perturbation theory may help [43| . The LO two-parton twist-3 contribution dominates in most range 
of Q^, and becomes comparable to the NLO twist-2 one as ~ 10 GeV^. The dominance comes from the chiral 
enhancement scale TOo(lGeV) = 1.74 GeV, which is of the same order of magnitude as low Q values. With this 
contribution being included, it is possible to accommodate the data in the large region. It supports the treatment 
in the PQCD approach to exclusive B meson decays [2H, in which the NLO twist-2 correction to B transition form 
factors is neglected, but the LO two-parton twist-3 one is included. That is, the latter should be the most important 
subleading contribution. The recent BaBar data for the pion transition form factor at low favors an asymptotic 
pion distribution amplitude [l^l . Figm-e [7] implies that the data of the pion form factor at low Q"^ also support 
the asymptotic models for the pion distribution amplitudes. Besides, the asymptotic models are favored from the 
viewpoint of light-cone sum rules, which reproduce the data of the pion form factor psj . As pointed out in p^ . the 
experimental data of F{Q^) at higher from [i^] have large additional uncertainties due to the model dependent 
substraction of transverse cross section, on top of the large statistical and systematic errors. Therefore, no conclusive 
statement on the agreement or disagreement between data and theory can be made. 

It is interesting to examine whether the shape of the pion distribution amplitudes matters for explaining the data, 
an issue that has been investigated and discussed in various aspects intensively. We test the following non-asymptotic 
pion distribution amplitudes |45| : 

<P^ix) ^ ^^x(l-.T)[l+0.16Cf^(u)+0.04Cr(z.)], 
^^{x) ^ ^-^[l + 0.59C,^/'(u) + 0.09Cy^«)], 

cl^lix) = -^[Cl/\u)+0.0l9Cl^\u)], (39) 



with the Gegenbauer polynomials 



Cy^u) = i(3u2 - 1), C^^iu) ^ ^{35u^ - 30^2 + 3), 

8 



1 

8 

15 



'(u) = §(5^2 - 1), Cl^\u) - y (21w^ - 14^2 + 1), (40) 

and the variable u = 1 — 2x. Figure [5] shows that the above models overshoot the data in the small region. Hence, 
the matrix elements of higher conformal spin operators may have been overestimated in (45j . 
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FIG. 9: Scale dependence of the pion form factor with ts{jJ.t) being defined in Eq. (|38|l . 



At last, we investigate the theoretical uncertainly in our calculation by studying the scale dependence, and the 
results are presented in Fig. [9l The curves corresponding to fif = t, fi = tg and fj.f = fi = t differ by about 10% in the 
considered range of Q^, implying that the theoretical uncertainty is under control. The sharp jump exhibited by the 
curve corresponding to fi{ — Q, fi = tg is attributed to the change of the QCD scale Aqcd with the upper bounds of 
the integration variables bi and 62. When the momentum transfer squared increases up to — 23 GeV^, the value 
A[^^q = 286 MeV is replaced by A^^^^^ = 204 MeV, and the sharp jump is produced as shown in Fig. [HI This jump 

is not obvious for /if = i (the other two curves), whose value is typically smaller than Q, so it does not excite Aq^q 

to at Q2 = 23 GeV^. 



IV. CONCLUSION 



In this paper we have extended the framework for higher-order calculations in the fc^ factorization theorem to 
the pion electromagnetic form factor. The prescription is that partons in both QCD quark diagrams and effective 
diagrams for TMD hadron wave functions are off mass shell by fcy. The light-cone divergences in naive definitions of 
TMD hadron wave functions were regularized by rotating the Wilson lines away from the light cone. Since the fc-p 
factorization theorem is appropriate for QCD processes dominated by contributions from small momentum fractions, 
terms suppressed by powers of momentum fractions have been neglected, and simple expressions can be obtained 
for the box and pentagon diagrams. We have demonstrated the disappearance of the soft logarithms, and the exact 
cancellation of the collinear logarithms in the difference of the quark and effective diagrams. Our work represents the 
first NLO analysis of the pion form factor in the fcy factorization theorem. 

Our numerical study has indicated that setting the renormalization and factorization scales to invariant masses 
of internal particles reduces NLO corrections. The NLO correction to the pion form factor becomes lower than 
40% of the LO one for > 7 GeV^. This observation supports the conventional scale choice adopted in the PQCD 
approach to exclusive processes based on the fc^ factorization theorem, and is consistent with that made in the collinear 
factorization theorem . We have found that the NLO twist-2 contribution does not help accommodating the data 
of the pion form factor. Instead, the LO two-parton twist-3 contribution plays a crucial role for the purpose. It is also 
confirmed that the asymptotic pion distribution amplitudes are favored by the data available at low . Comparing 
the results from the different choices of the renormalization and factorization scales, the theoretical uncertainty in 
our calculation is estimated to be around 10%. We shall apply the same framework to the NLO analysis of B meson 
transition form factors in a forthcoming paper. 
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